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Abstract 

The presence of losses in nonlinear photonic structures is a crucial issue for modern applications. 
Active parts are introduced for wave power compensation resulting in unbalanced gain and loss 
landscapes where localized beam propagation is, in general, dynamically unstable. Here we provide 
generic sufficient conditions for the relation between the gain-loss and the refractive index profiles 
in order to ensure efficient wave trapping and stable propagation for a wide range of beam launching 
conditions such as initial power, angle of incidence and position. The stability is a consequence of 
an underlying dynamic power balance mechanism related to a conserved quantity of wave dynamics. 
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Modern photonic applications utilize the combination of nonlinearity and inhomogeneity 
in order to provide advanced functionality in properly engineered metamaterials and metade- 


vices 




2| These structures have the form of multilayered media consisted of materials such 


as ordinary dielectrics as well as metals and graphene. The presence of metals results in 
plasmonic excitations that can boost nonlinear effects due to high field values and small 


mode volume 


wherea. the pteeence of gtaphene layete ie aeeompan.ed whh ve^ ettong 


Kerr nonlinearities, jd] both resulting in the formation of self-localized modes. |^l-lsl] The 
nonlinearity plays a crucial role in functionality related to dynamic and all optical light con¬ 
trol through wave-material and wave-wave interactions. However, both ordinary dielectrics 
and metal or graphene layers introduce significant losses that can hamper the nonlinear 


i. B 


functionality of these structures by restricting the wave propagation to small distances. 
This crucial drawback of the respective photonic structures necessitates the utilization of 
active parts (hot-spots) in the form of doped and pumped dielectrics in order to provide 
the necessary gain for loss compensation, liJQ introducing an inhomogeneous gain-loss 
landscape. A similar type of nonconservative inhomogeneity also appears in applications 


related to soliton-forming laser cavities 


ities. jl^ 


In all these cases the formation and robust 


propagation of a self-localized mode is determined by both the diffraction-nonlinearity and 
the loss-gain balance, which cannot be considered separately. 

From an engineering aspect, even the presence of spatially homogeneous gain and loss 
in an optical lattice can significantly enrich soliton dynamics providing soliton routing and 
acceleration functionalities, jl4, ^ in addition to the numerous applications related to con¬ 
servative lattices that include the formation of solitons, surface waves and defect modes. 
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lilt Moreover, the appropriate design of gain and loss inhomogeneity provides another 


degree of freedom for wave manipulation. 




The symmetry properties of the inho¬ 


mogeneity profiles have been shown to play a crucial role on the system features. It has 
been shown that, for the case of VT symmetry, the system has a real spectrum and sup¬ 
ports a continuous family of solitons. (20|, l2lt Other types of symmetries that restrict, not 
the profiles of the refractive index and the gain-loss inhomogeneity as in the VT case, but 
i;heir mutual relation, have been also shown to support such continuous soliton families, 
0, 22, 22| in contrast to the common case of dissipative solitons where, in general, only 
isolated solitons exist. 

A solitary wave can propagate at a fixed transverse position of a planar structure near 
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the interface between a gain- and a loss-region, where a static power balance condition is 
satisfied. However, any deviation from this specific position or from a zero angle of incidence 
can lead to continuous power increasing or decreasing, resulting in an unstable behavior, as 


in the case of stationary solitons pinned to hot-spots. 


24142811 The utilization o 


modulations of the linear or the nonlinear refractive index has been proposed 
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spatial 
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introducing effective potential wells resulting to wave trapping in the specific position and 
preventing large excursions within the two regions around the fixed position. Even in such 
cases, wave oscillations around the balance position, can be unstable when the gain and loss 
of the interfaced parts are unbalanced, as in the most typical case where narrow hot-spots 
with high gain are utilized in order to compensate for more extended parts with relatively 
small losses. The instability arises from the fact that the dynamic power balance of the 
wave depends on the extent of the oscillations in the two parts, since the wave amplification 
and attenuation in the two phases of the oscillation are not equal in general. Therefore, an 
appropriate refractive index modulation has to take into account the gain and loss profile, 
in order to ensure a dynamic balance of gain and loss and a stable wave propagation. 

In the following, we present generic efficient conditions for the relation between the gain- 
loss and the refractive index profiles allowing, not only for stable stationary propagation 
at a specific point with a zero angle of incidence, but also for dynamic power balance for 
localized waves with a wide range of positions and angles of incidence, that are applicable to 
any type of planar photonic structure that may have unbalanced gain and loss properties. 


Model and Method 

Nonlinear wave propagation in a transversely inhomogeneous planar photonic structure 
is described by the NonLinear Schrodinger Equation (NLSE) 

iuz + Uxx + \y{x) — iW (x)] u + 2|upu 0 (1) 

where u is the normalized electric field envelope, and x are the normalized longitudinal 
and transverse dimensions, and V{x), W{x) are the transverse refractive index and gain- 
loss profiles, respectively. Eor spatially localized (solitary) waves, we can define the useful 
quantities corresponding to the wave mass m = f \u\^dx and momentum p = i f (uu* — 
UxU*)dx. In the absence of inhomogeneity (V = W = 0), m and p are conserved, whereas 
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under the presence of inhomogeneity they vary as 

dm 

dz 
dv 


m— = 
dz 


= r(xo) 


dXt] 


where 


r(xo) = 2 


r*+oo 


\u{x — Xq)'\^W {x)dx 


^effi^o) — 2 


r‘ + 00 


\u{x — Xo)\‘^V {x)dx 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


are the mass variation rate and the effective potential, respectively, and xq is the wave center 
varying as dxo/dz = p/m = v with the velocity v corresponding to the propagation angle. 
Therefore, the solitary wave propagates as an effective particle with mass and momentnm 
variations depending on the nonconservative [lT(x)] and the conservative [l^(x)] part of 
the inhomogeneity, respectively. Wave propagation dynamics are described in the three- 
dimensional space {xo,v,m). It can be readily shown 

dV{x 


19| that, nnder the condition 


dx 


- C'IT(x), 


( 6 ) 


the existence of an exact invariant of the motion, given hj K = Clnm + v, restricts the wave 
dynamics in a two-dimensional snrface. This is a general property of any type of solitary 
wave in the presence of inhomogeneities of arbitrary profile and magnitude. The condition 
([HD ensures the static power balance for a stationary solitary wave located at a fixed point 
r(xo) = 0 at the vicinity of the interface between a lossy and an amplifying part. Moreover, 
it is a stronger condition, sufficient for the dynamic power balance of solitary waves with 
nonzero angles of incidence and positions deviating from the fixed point that undergo stable 
oscillations, as we show in the following. Notice that the condition (ED is qnalitatively 
different from the VT symmetry condition, since it does not imply any restriction on the 
symmetry properties of the nonconservative [lT(x)] and the conservative [V(x)] part of the 
inhomogeneity, bnt only a mutual relation of their profiles; therefore, it is also applied in 
non-symmetric profiles. Under condition (pD when V{x) is even, W{x) is odd (and vice 
versa) as in VT symmetric confignrations. However, the VT symmetry condition suggests 
only the existence of a fixed point and does not ensure its stability. 

In the following, we exploit the consequences of this condition with respect to the dynamic 
power balance for solitary waves in a wide variety of planar structures, and prove that Eq. (jED 
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serves as a generic sufficient condition for stable wave propagation in accordingly designed 
photonic structures. 


Results and Discussion 


We focus on multilayer photonic structures with piecewise constant gain and loss prohles. 
According to Eq-lEI), the linear refractive index profile is a piecewise linear function, so that 


W{x) = V{x) = A*(^) 

i i 

with 


( 7 ) 


ni(a:) 


Ai{x) 


^i,l ^ X ^ Xi^2 

0, elsewhere 

CiX + di, Xi^i < X < Xi^2 

0, elsewhere 


( 8 ) 

( 9 ) 


The dynamics of solitary wave propagation in such structures is determined by Eqs. 
(El), (El)- Without loss of generality, in order to simplify our analysis and provide intuitive 
understanding, we consider inhomogeneities of relatively small amplitude so that we can 
obtain closed form equations for the Eqs. (Ill)-(l5]) by utilizing the soliton solution of the 
homogeneous NLSE u = risech[ri{x — xq)] exp[i(n/2)x + i{rf‘ — n^/4)^] in the calculation of 
the respective integrals, resulting in 


r(a:o) = m [TTj (xi,2) - TTj (xi,i)] 

i 

U^ffixo) = -2m^[Xi{xi^2) - 

i 


( 10 ) 

( 11 ) 


7ri(x) 

Xi{x) 


tti tanh —ix 
2 ^ 

Ci{x - Xq) - 


- Xo) 

CiX + di 
^m{x-xo) 2 


— In + 1) 

m 


( 12 ) 

( 13 ) 


with m = 2r]. 

Eirst, we consider a structure consisting of two interfaced semi-infinite parts with unequal 
gain and loss coeffecients and linear refractive indices profiles fulfilling the condition dH]) with 
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C* = — 1 as shown in Fig. 1(a). The effective potential has a local minimum, corresponding 
to a fixed point, in the vicinity of the interface with its exact poisition depending on the 
soliton mass. For a soliton of mass m = 1 the fixed point is located at Xq = —0.69. Stable 
propagation of a stationary soliton with initial position at the fixed point is shown in Fig. 
1 (b). The fulfillment of the condition (Ej) results in refractive index slopes appropriately de¬ 
fined in terms of the gain and loss coefficients in each part. In terms of soliton dynamics, the 
direct consequence of the condition is that the trapping potential is such that no continuous 
mass increase or decrease takes place as the travelling distance in each region is such that 
the soliton spends less time in the high gain region than in the low loss region. In fact, this 
dynamic power balance mechanism results in asymptotic evolution to the stable fixed point 
(attractor). As shown in Fig. 1(c), the effective particle orbit for a soliton initially located at 
Xo = —10 evolves to the fixed point, while remaining in the aforementioned two-dimensional 
surface of the phase space. The rate of convergence to the fixed point orbit increases with 
the magnitude difference between the gain and loss coefficients and for the specific case is 
quite small as shown in Fig. 1(d). Notice that the period of oscillations scales with 
The importance of the condition (ED can be shown in comparison to the cases where it is 
not fulfilled, resulting in either continuous mass decreasing or increasing as shown in Fig. 
1(e) and (f), respectively, and unbounded phase space orbits [Fig. 1(c)]. 

A typical realistic case with practical importance is a planar structure consisted of an 
amplifying part of finite width (hot-spot) in a lossy medium, with gain-loss and refractive 
index profiles fulfilling the condition (ED as shown in Fig. 2(a). For a soliton of mass m = 1 
and the parameters values of the specific structure, an asymmetric potential well can be 
formed as shown in Fig. 2(a). The fixed point located at Xq = —6.25 corresponds to 
stable stationary soliton propagation [Fig. 2(b)] whereas stable large amplitude oscillations 
(asymptotically evolving to the stationary soliton) can take place as a consequence of the 
condition for dynamic power balance, as shown in Fig. 2(c). Since the potential well has 
not an infinite depth, initial soliton conditions corresponding to untrapped dynamics result 
in travelling solitons of continuously decreasing mass, as illustrated in Fig. 2(d). It is worth 
noticing that the trapping conditions depend on both the parameters of the structure and 
the soliton mass, so that in each structure solitons having a mass below a critical value 
cannot be trapped. This fact results from the interplay of the two spatial scales, namely the 
soliton width (~ m~^) and the amplifying part width (Ax) as well as the relative magnitude 
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of the gain and loss coefficients and is reflected in the effective particle model as a bifurcation 
of the fixed point corresponding to the local minimum of the effective potential. The spatial 
width and the depth of the potential well, for a given soliton mass value, determine the range 
of initial positions and angles of incidence (velocities) for soliton trapping and stability. 

The case of a structure with two hot-spots is considered in Fig. 3. As expected, the 
increased complexity of the structure results in richer dynamics and trapping capabilities. 
In such case, we can have two fixed points, as shown Fig. 3(a). Therefore, under dynamic 
power balance conditions, trapping and stable soliton oscillations can occur either on the left 
potential well [Fig. 3(b)] or on the right one [Fig. 3(c)]. Moreover, for appropriate initial 
conditions, extended stable oscillations can take place in the region above the two potential 
wells, for effective particle energy above the left and below the right local maximum of the 
effective potential [Fig. 3(d)]. The existence and bifurcations of the two fixed points again 
depend on the soliton mass and the parameters of the structure, so that we can have two, 
one or zero fixed points for a given soliton mass value. 

Robust coexistence of two solitons trapped in the two different potential wells is shown 
in Fig. 4(a) for soliton mass m = 1. Although initially located at the corresponding fixed 
points, solitons oscillate due to mutual interactions depending strongly both on the soliton 
width (mass) and the distance between the two hot-spots. Different interaction scenarios are 
possible as in the case of two solitons of higher mass {m = 1.5) as shown in Fig. 4(b) where, 
although both solitons are launched at the corresponding fixed points, under interaction 
the right one is detrapped and travels with continuously decreasing mass whereas the other 
is stably trapped in its potential well. A numerous list for interaction scenarios, including 
solitons of different masses, can be considered which can be very interesting in terms of 
light control applications. It is worth emphasizing that it is the fulfillment of the dynamic 
power balance condition between the refractive index and the gain-loss profiles that allows 
for stable soliton dynamics and mutual interactions that could not take place either in the 
absence or in the inappropriateness of the trapping potential. 

Conclusions 

The fundamental problem of power balance of a nonlinear wave in a photonic structure 
with unbalanced gain and loss has been addressed. Sufficient conditions between the re- 
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fractive index and gain-loss profiles have been derived for dynamic power balance of soliton 
propagation. In contrast to static power balance, that ensures only the existence of a fixed 
point corresponding to stationary soliton propagation, the dynamic power balance ensures 
the stability of the fixed point solution, allowing for stable propagation for a wide range of 
initial soliton positions and velocities, which is crucial for realistic applications. The anal¬ 
ysis has been based on a simple effective particle model providing, not only the sufficient 
conditions, but also intuitive understanding of the complex soliton dynamics and being in 
remarkable agreement with the full model. The concepts and results of the dynamic power 
balance, illustrated here for simplicity only for piecewise constant gain-loss profiles, are so 
general that can be directly applied to any type of gain-loss profiles, nonlinear refractive in¬ 
dex and nonlinear gain-loss inhomogeneities, as well as two-dimensional photonic structures. 
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(e) (f) 

FIG. 1: Single interface between a lossy (ai = —0.0005) and an amplifying (a 2 = 0.001) region, 
(a) Gain-loss W{x) and refractive index V{x) profiles [top]; Effective potential Ueff{xo) and mass 
variation rate T(xo) for a soliton of mass m = 1 for a refractive index profile fulfilling the condition 
(|6j) with C = —1 {ci = —ai^ i = 1, 2) [bottom], (b) Stationary propagation for initial soliton position 
xq = —0.69 corresponding to the fixed point depicted by a thick dot in (a), (c) Phase space orbits 
of a soliton with initial position xq = —10 under dynamic balance conditions q = —= 1,2 
(blue), and for unbalanced cases with ci = —ai,C 2 = —7a2 (red), ci = —ai,C 2 = —0.25a2 (green). 
(d),(e),(f) Soliton propagation for conditions corresponding to the three orbits shown in (c). The 
thick black lines depict results from the effective particle model. 
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(C) (d) 

FIG. 2: Single hot-spot structure consisted of an amplifying part (a 2 = 0.001) of finite width 
Ax = 4 in a lossy medium (ai = as = —0.0005) for a refractive index profile fulfilling the condition 
with C = —1 {ci = —ai^i = 1,2). (a) Gain-loss W{x) and refractive index V{x) profiles [top]; 
Effective potential t/e//(xo) and mass variation rate r(xo) for a soliton of mass m = 1 [bottom], 
(b) Stationary propagation for initial soliton position xq = —2.65 corresponding to the fixed point 
depicted by a thick dot in (a), (c) Trapped soliton oscillations for initial position xq = 1. (d) 
Travelling soliton propagation for initial position xq — —10. The thick black lines depict results 
from the effective particle model. 
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(C) (d) 

FIG. 3: Double hot-spot structure consisted of two amplifying parts (a 2 = — 0.001) of finite 

width Ax = 4 located at Xc = ±5 in a lossy medium (ai = as = — 0.0005) for a refractive 

index profile fulfilling the condition (|6j) with C = —1 (q = —= 1,2). (a) Gain-loss W{x) and 
refractive index V{x) profiles [top]; Effective potential t/e//(xo) and mass variation rate T{xo) for a 
soliton of mass m — 1 [bottom], (b) Trapped soliton oscillations in the left potential well for initial 
position xq = —11. (c) Trapped soliton oscillations in the right potential well for initial position 
xq = 4. (d) Extended trapped soliton oscillations for initial position xq = —14. The thick black 
lines depict results from the effective particle model. 
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FIG. 4: (a) Oscillatory interaction of two solitons of mass m — 1 and initial positions xq = —7.65 
and xo = 2.33 corresponding to the two fixed points depicted in Fig. 3(a). (b) Interaction of two 
solitons of mass m = 1.5 and initial positions xq = —7.46 and xq = 2.54 corresponding to the 
respective fixed points; the left soliton remains trapped whereas the right soliton is detrapped. 
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